arXiv:1503.05895v3 [astro-ph.CO] 7 Aug 2015 


Cosmic microwave background radiation temperature in a dissipative universe 

Nobuyoshi Komatsu^B and Shigeo KiMURA^ 

^Department of Meehanical Systems Engineering, Kanazawa University, 

Kakuma-machi, Kanazawa, Ishikawa 920-1192, Japan 
^ The Institute of Nature and Environmental Technology, 

Kanazawa University, Kakuma-machi, Kanazawa, Ishikawa 920-1192, Japan 
(Dated: August 10, 2015) 

The relationship between the cosmic microwave background radiation temperature and the red- 
shift, i.e., the T-z relation, is examined in a phenomenological dissipative model. The model contains 
two constant terms, as if a nonzero cosmological constant A and a dissipative process are operative 
in a homogeneous, isotropic, and spatially flat universe. The T-z relation is derived from a general 
radiative temperature law, as appropriate for describing nonequilibrium states in a creation of cold 
dark matter (CCDM) model. Using this relation, the radiation temperature in the late universe is 
calculated as a function of a dissipation rate ranging from /I = 0, corresponding to a nondissipative 
ACDM model, to fi = 1, corresponding to a fully dissipative CCDM model. The T-z relation for 
/t = 0 is linear for standard cosmology and is consistent with observations. However, with increas¬ 
ing dissipation rate fi, the radiation temperature gradually deviates from a linear law because the 
effective equation-of-state parameter varies with time. When the background evolution of the uni¬ 
verse agrees with a fine-tuned pure ACDM model, the T-z relation for low jl matches observations, 
whereas the T-z relation for high fi does not. Previous work also found that a weakly dissipative 
model accords with measurements of a growth rate for clustering related to structure formations. 

These results imply that low dissipation is likely for the universe. The weakly dissipative model 
should be further constrained by recent observations. 

PACS numbers: 98.80.-k, 98.80.Es, 95.30.Tg 


I. INTRODUCTION 

The blackbody radiation temperature of the cosmic 
microwave background is Tq = 2.725 ± 0.002 K at the 
present time [I| . This temperature is regarded as the ev¬ 
idence for a hot Big Bang. That hypothesis is supported 
by measurements of the radiation temperature-redshift 
T-z relation nn. The observations are consistent with 
a linear law, T = To(l -|- z), as can be derived from stan¬ 
dard cosmology, such as ACDM (lambda cold dark mat¬ 
ter) models. However, the ACDM model has several the¬ 
oretical difficulties [l4 - [l^ . although it can elegantly ex¬ 
plain the accelerated expansion of the late universe [H- 
[23|. To explain the acceleration, various models have 
been suggested, such as A(t)CDM which assumes a time- 
varying cosmological term [23 - 1^ . bulk viscous models 
which assume a bulk viscosity for the cosmological fluid 
[29l - l3^ . and CCDM which assumes the creation of cold 
dark matter 

The bulk viscous and CCDM models presume the ex¬ 
istence of irreversible entropy in a homogeneous isotropic 
universe, unlike the ACDM and A(t)CDM models. For 
example, in the CCDM model, irreversible entropy is gen¬ 
erated from gravitationally induced particle creation in 
nonequilibrium thermodynamic states [13, H^. A pos¬ 
sible equivalence of the bulk viscosity and matter cre¬ 
ation dissipative mechanisms has been discussed in Ref. 
[ 3 ^ . In addition, the connections between warm infla- 
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tion [ 4 ^, A(t)CDM models, and CCDM models have 
been debated in Ref. [s^- To examine the dissipative 
processes, a general radiative temperature law for adia- 
batic particle creation has been proposed by Lima et al. 
[33l - l36j| . That law has been examined from various view¬ 
points [36l - [38l| . In particular, a simple T-z relation of 
the form T = To(l -|- z)^~^ is frequently compared with 
observations, where ,9 is a constant parameter (cf. Refs. 
0111). The simple T-z relation is obtained from a gen¬ 
eral radiative temperature law, if the effective equation- 
of-state parameter We is constant [H, HI. However, in 
the CCDM model, We varies during the evolution of the 
universe. Thus, a time-varying We needs to be consid¬ 
ered when the radiation temperature is discussed in the 
CCDM model. However, the radiation temperature in 
the CCDM model has not yet been quantitatively exam¬ 
ined from this viewpoint. It is important to do so in order 
to acquire a deeper understanding of the CCDM model. 

In the CCDM model, a negative sound speed [11 and 
the existence of clustered matter [dl are necessary to 
properly describe the growth rate for clustering related 
to structure formations. Alternatively, a phenomenolog¬ 
ical dissipative model M has been proposed, in which 
the entropic force [Hsl is modified. The model as¬ 
sumes constant terms that are equivalent to a nonzero 
cosmological constant A and a dissipative process. In 
previous work 0|, the dissipation rate was varied from 
jJ — D, corresponding to a nondissipative ACDM model, 
to ft = 1, corresponding to a fully dissipative CCDM 
model. Low dissipation was found to correctly describe 
observations of structure formation. The dissipation rate 
pL is expected to affect the T-z relation, not only because 
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the effective equation-of-state parameter We depends on 
/i, but also because We varies during the evolution of the 
universe. 

This model makes it possible to examine a dissipa¬ 
tive universe systematically, ranging from a nondissipa- 
tive ACDM model to a fully dissipative CCDM model. 
To clarify the properties of the radiation temperature in 
a dissipative universe, the T—z relation can be examined. 
A general radiative temperature law is applied to the dis¬ 
sipative model, to formulate the T-z relation. Based on 
this formulation, the radiation temperature in the late 
universe is calculated numerically as a function of the 
rate of dissipation. The present study provides new in¬ 
sights and a unique approach for examining a dissipative 
universe. 

The remainder of the article is organized as follows. 
In Sec. im the general radiative temperature law for adi¬ 
abatic particle creation is briefly reviewed. In Sec. Ell 
a phenomenological modified dissipative model is pro¬ 
posed. In Sec. IIYI the temperature law is applied to 
the modified model, and the T-z relation is formulated. 
In Sec. El the radiation temperature in a dissipative uni¬ 
verse is examined. Finally, in Sec. ED the conclusions are 
presented. (The modified dissipative model is different 
from dissipative particle dark matter models examined 
in Ref. [s^.l 

II. GENERAL RADIATIVE TEMPERATURE 
LAW FOR ADIABATIC PARTICLE CREATION 


where n, s, e, and p are the particle number density, en¬ 
tropy density, energy density, and pressure, respectively. 
Here, T and Pc are the particle production rate and the 
dynamic creation pressure, respectively. When T = 0 and 
Pc = 0, the three balance equations reduce to the conser¬ 
vation law for equilibrium states in a standard cosmology. 
Equation Q can be rewritten as N/N = T, using the 
total number N oc na^ of particles in the comoving vol¬ 
ume [ 3 ^ . Keep in mind that the entropy per particle 
a = S/N is assumed to be constant [H, [s^, i.e., ct = 0, 
where S oc sa^ is the entropy in the comoving volume. 
The constant value of a indicates SjS = N jN = T, 
which has been used for calculating the right-hand side 
of Eq. dH). (Units are chosen so that c = fcs = I, unless 
otherwise stated. Here ks is the Boltzmann constant.) 

The local Gibbs relation should be valid even for the 
nonequilibrium process considered here [s^. Accord¬ 
ingly, the thermodynamic quantities are related to the 
temperature T by 

nksTd (—^ = nksTda = de --^dn. (6) 

Vn/ n 

Substituting <7 = 0 into Eq. using both this result 
and Eq. ®, and rearranging Eq. one obtains the 
dynamic creation pressure 

Pc = -{e+p)-^. (7) 

The thermal evolution of matter creation can then be 
written as 


In this section, the general radiative temperature law 
for adiabatic particle creation is reviewed, following the 
work of Lima et al. A homogeneous, isotropic, 

and spatially flat universe is initially considered. The 
line element given by the Friedmann-Robertson-Walker 
(FRW) metric is 


ds^ = c^df - a'^{t){dr'^ + r'^dO'^ + sin^ 9d4>^) (1) 


where c is the speed of light and a{t) is the scale factor 
at time t. The Hubble parameter H is 



( 2 ) 


In what follows, nonequilibrium thermodynamic states 
of the cosmological fluid in a FRW background are con¬ 
sidered. To this end, adiabatic particle creation is as¬ 
sumed The balance equations for the number of 

particles, entropy, and energy [36| can then be written as 


n + SHn = nF, 


(3) 


s -I- 3Hs = sF, 


(4) 


and 


T / dp's h 

T ~ \de)^ n 



(F - 3H). 


( 8 ) 


This equation has been derived by Lima et al. 

(More general formulations have been examined in Ref. 
[33| . Harko has recently developed an equivalent formula¬ 
tion for a modified gravity theory with geometry-matter 
coupling [131.) 

For the nonequilibrium thermodynamic states consid¬ 
ered here, a radiation temperature relation can be ob¬ 
tained from Eq. ([8]). Substituting p = ejS into Eq. ([8]), 
one finds 


T a r /r j- • ^ 

— =-1— (for radiation). 

T a 3 


That can be rearranged as 
1 d{aT) F 


aT dt 


diaT) F , 
- or i = -dt. 

3 aT 3 


(9) 


( 10 ) 


Integrating Eq. m from arbitrary time t to the present 
time to gives 


r°'^° d{a'T') 

laT 


T{t'] 


'T' 


dt' 


( 11 ) 


and therefore 




e -I- 3H{e -I- p + pc) = 0, 


(5) 


In 


( 12 ) 
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where oq and Tq are the present values of the scale fac¬ 
tor and of the radiation temperature, respectively. This 
equation can be rearranged as 


T = To exp 



(13) 


so that 


T = Tn 


or equivalently 




T = To(l -I- z) exp 


1 rv! ^' 


(14) 


(15) 


where the normalized scale factor a and the redshift 2 
are 


a ao . 

a = — and z =-1. 

ao a 


(16) 


Equations (nsi), (m, and (USD are the general radia¬ 
tive temperature law for adiabatic particle creation [3^ . 
In this paper, that law is used to calculate the radi¬ 
ation temperature T in a dissipative universe. When 
r = 0, Eq. (IT5|) reduces to the linear law T = ro(l -|- z). 
However, T is not zero for the nonequilibrium states 
considered here. Therefore, it is necessary to compute 

T{d'){dt'/da')da' in Eq. (fT4l) and T{z'){dt'/dz')dz' 
in Eq. m to calculate T. These integrals depend on the 
limits, the particle production rate, and the background 
evolution of the universe. In Sec. El the calculated T 
is compared with observations in the late universe. Ac¬ 
cordingly, in the present study, the limits of integration 
in Eqs. (II31)-(1I5D correspond to the late universe. 

A simple temperature relation is obtained from Eq. ® 
for specific cases [H, . For example, if T = 3,5i4 = 

3/3(d/a), integrating Eq. ([S]) gives 


r = To(^)' = To(l + z)i-^ (17) 

where P is taken to be constant. This simple temperature 
relation has been studied in detail [s^. In contrast, here 
a general radiative temperature law is considered for a 
dissipative universe. 


III. MODIFIED DISSIPATIVE MODEL WITH 
CONSTANT TERMS 

Entropic cosmology has been proposed to explain 
the accelerated expansion of the universe [s^. In the 
entropic-force model, the horizon of the universe is as¬ 
sumed to have an associated entropy and an approximate 
temperature due to information holographically stored 
there [s^. Recently, various entropic-force models have 
been examined in detail IblMbSlI . For example, in place 
of the Bekenstein entropy the Tsallis-Cirto entropy 


61j based on nonextensive statistics l62l| h as been ap¬ 
plied to the horizon of the universe [571 l58l |. Basilakos 
et al. ( 53 . [ 5 ^ have shown that simple combinations 
of pure Hubble terms, such as H^, H, and H, are in¬ 
sufficient for a complete description of the cosmological 
data (5ll |. Therefore, the constant term plays an im¬ 
portant role. Accordingly, a phenomenological entropic- 
force model that includes constant terms, an irreversible 
entropy ^irr, and a kind of reversible entropy Srev has 
been proposed [5l|. Using this model, a dissipative uni¬ 
verse can be analyzed systematically, over the entire 
range from a nondissipative ACDM model to a fully dissi¬ 
pative CCDM model. In the present work, it is called the 
modified dissipative model. In what follows, it is briefly 
reviewed according to Ref. [5l|. Keep in mind that the 
entropic-force considered here is different from the idea 
that gravity itself is an entropic force [^, . 

The Friedmann, acceleration, and continuity equations 
for the modified dissipative model become 


2 SttG 2 

H = -—p + aHl, 


(18) 


d 

a 


47rG 

47rG 


(^P+ + lirrHo 

(^p+ -h aiLp, 


and 



with two dimensionless constants 

a > 0 and qirr > 0. 


(19) 


( 20 ) 

( 21 ) 


Here G, Hq, and p are the gravitational constant, the 
Hubble parameter at the present time, and the mass den¬ 
sity of the cosmological fluid, respectively [5l|. The mass 
density is p = ejt?. The effective pressure Pe in Eqs. m 
and (1^ is 

Pe = P + Pc (22) 


where pc is a pressure derived from irreversible entropy 
related to dissipative processes. In this study, Pc is taken 
to be equivalent to the dynamic creation pressure in the 
CCDM model. (For A(t)CDM models, nonzero terms 
related to reversible entropy appear on the right-hand 
side of the continuity equation . That model involves 
the transfer of energy between two fluids [6^ - 1^ . 1 

The modified dissipative model is used to calculate a 
radiation temperature T from the general radiative tem¬ 
perature law in Eqs. (nnii-disi)- In Sec. El the T-z re¬ 
lation in the late universe is examined, to compare the 
calculated value of T with observations. The pressure 
of the cosmological fluid in the present model is negli¬ 
gible, p = 0. A matter-dominated universe is assumed. 
Consequently, the effective pressure Pe is 


Pe = p -I- Pc = Pc = - 


47rG 


7irr- 


(23) 
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On the other hand, substituting pc = Pe, P = 0, and 
e = pc? into Eq. one obtains 


where is a constant parameter related to dissipative 
processes, 



Accordingly, from Eqs. (l23)) and (l24ll . 


3g girr^o 

47rG p 


(24) 


(25) 


Assume that the oiH^ terms in Eqs. m and (HU) are 
equivalent to the cosmological constant A/3 in the stan¬ 
dard ACDM model. That is, dSr^y = 0 for homogeneous 
systems in CCDM models [^|^. In contrast, Pe is re¬ 
lated to the irreversible entropy in dissipative processes 
because Pe = Pc- In other words, the term in 

Eq. (fT^ is related to the irreversible entropy. Accord¬ 
ingly, the preceding cosmological equations are equiva¬ 
lent to those for an extended ACDM model in a dissi¬ 
pative universe. Therefore, the ai?Q term is interpreted 
as a modification of the Einstein tensor. In contrast, Pe 
is a modification of the energy-momentum tensor of the 
Einstein equation M- The properties of this model are 
equivalent to those of the modified entropic-force model 
examined in Ref. M- 

Combining Eqs. m and (HU), and using Eq. (l23l) . one 
finds 


H = -CmH^ + CagHl (26) 

where Cm and Cag are dimensionless constants [5l| 

Cm = l.h and Cag = ^^^^. (27) 

Here Cm = 1.5 corresponds to a matter-dominated uni¬ 
verse in standard cosmology EH- Solving Eq. (l26t . 
one obtains 

= (l-^2A)d-3 + DA 

= tlmtt ^ -j- Da (28) 


where d = a/ao in Eq. (1161) . The two constant parameters 
are defined as 

Da = and D^ = 1 — Da. (29) 


This solution is the same as that in the standard ACDM 
model [m, . Accordingly, the constant term Da be¬ 

haves as if it were Da in the standard ACDM model. 
Similarly, D^ behaves as if it were D^- Here D^ and 
Da are the density parameters for matter and for A, re¬ 
spectively. The density parameter D,. for radiation is 
neglected in the late universe. 

To study a dis^ative universe quantitatively, define a 
dissipation rate |51 | 


7irr _ Dg 
Cag Da 


(30) 


Od 3 (31) 

As discussed in Ref. [5l|, when yin- = 0, one obtains 
/t = 0 from Eq. ((30ll . In this case, the present model is 
equivalent to the standard nondissipative ACDM model. 
In contrast, when a = 0, one obtains p = 1 from Eq. dSQl) 
because Cag = = yjj.,.. in that case, the present 

model is equivalent to the fully dissipative CCDM model 
proposed by Lima et al. [4l|. That is, p = 0 corresponds 
to the nondissipative ACDM model, whereas p = 1 cor¬ 
responds to the fully dissipative CCDM model. In this 
way, the extent of the dissipative universe is determined 
by the dissipation rate p. 

Next, consider the effective equation-of-state parame¬ 
ter in the modified dissipative model. From Eq. ()24[) . We 
becomes 


We = 



r 

'm' 


After some algebra, one finds 


We = 


Cljgd^ 

1 — Da + Ddu^ 


Dm + Dijd^ 


(32) 


(33) 


where d is the normalized scale factor a/ag. For details, 
see Ref. [5l|. Note that We is not equal to the equation- 
of-state parameter w for a generic component of matter, 
because w is zero in a matter-dominated universe, for 
which p = 0. (A more general inhomogeneous equation 
of state has been examined in Ref. [^.1 


IV. RADIATION TEMPERATURE IN THE 
MODIFIED DISSIPATIVE MODEL 


The general radiative temperature law in Eq. (HI is 
rearranged as 


T = To exp 
= To(l -I- z) exp 


df' 

T{d')—dd' 
^ aa' 


4 /'r(d')^rfa' 

3 A da' 


(34) 


First consider the calculation of dt/dd in Eq. (1341) . (For 
simplicity, the prime is omitted.) As discussed in Sec. 
nn the background evolution in the modified dissipative 
model is the same as that in the standard ACDM model. 
Accordingly, solutions of the standard ACDM model are 
used for a spatially flat universe. The solution is 


d{t) 


/ - \ 1/3 

( ] sinh2/3 

\nA J 


3i7( 




(35) 
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V. EVOLUTION OF THE RADIATION 
TEMPERATURE IN A DISSIPATIVE UNIVERSE 


t{a) = 


SHoVfh 


sinh 


-1 



(36) 


The density parameter for radiation is negligible in this 
late universe. Differentiating Eq. (1551) with respect to 
t, substituting Eq. (1551) into the result, and rearranging, 
one obtains 


dt sinh^'^ q(d) cosh ^ g(d) 


where 



(37) 


(38) 


Next, rearrange E in Eq. (1551) . Using Eq. (1551) . it can be 
written as 

r = -3HWe = -3 HoWe. (39) 

Substituting Eqs. (l28l) and (1551) into Eq. (1591) leads to 


r = 


3i7o^rf\/Umd ^ + Ua 


(40) 


Substituting Eqs. dSil) and (l40l) into Eq. (1541) , one finds 


T =To(l + z) exp 
=ro(l + z) exp 

where 


1 




3 

—D d 




K{d')dd' 


(41) 


= ^^'^ -sinh^/^g(d)cosh ^g(d) (42) 

a LjYid “1“ i 2 


and g{d) is given by Eq. (1551) . Equation (15T1) is the radi¬ 
ation temperature-redshift relation for the modified dis¬ 
sipative model in the late universe. The infiuence of dis¬ 
sipation is included in Using Eq. m, the radiation 
temperature can be determined as a function of the dis¬ 
sipation rate g, where g — Ud/Da from Eq. (1501) . 

As a specific case, a simple temperature relation can 
be obtained from Eq. For example, if We = —T/ {3H) 
from Eq. (1321) is assumed to be constant, then 

T = Toil + zY^^^. (43) 


This equation is equivalent to Eq. dm, replacing We by 
—j5. In general, We is not constant in a dissipative uni¬ 
verse, as examined in the next section. Accordingly, Eq. 
m plays an important role in studying the dissipative 
universe. 


The dissipation rate p, from Eq. dsni) is 


7irr _ Hd 

Cag Da 


(44) 


Zero dissipation corresponds to a nondissipative ACDM 
model, whereas p = 1 corresponds to a fully dissipa¬ 
tive CCDM model. As discussed in Ref. [^, Da can 
be determined from the background evolution of the uni¬ 
verse. Accordingly, Da = Da from a fine-tuned standard 
ACDM model. Consider a spatially fiat universe in which 
(Dm, Da) = (0.315,0.685) based on the Planck 2013 re¬ 
sults [13. That is, let Da = Da = 0.685 [5l|. To confirm 
the background evolution of the universe in the present 
dissipative model, consider the luminosity distance dz, 
[nl given by 




dy 

FivY 


(45) 


The integrating variable y and the function Fiy) are 
given by 

y = — = d~^ and Fiy) = ^ (46) 

a Ho 


1.0 
0.5 

^ 0.0 

^ -0.5 

f -1-0 

-1.5 

- 2.0 

0.0 0.5 1.0 1.5 2.0 

Z 

FIG. 1: (Color online). Dependence of the luminosity dis¬ 
tance dL on the redshift «. The three continuous curves rep¬ 
resent the modified dissipative model for Da = 0, 0.685, and 
1. They are respectively equivalent to the standard ACDM 
model for (Dm, Da) = (1, 0), (0.315, 0.685), and (0,1) in a spa¬ 
tially fiat universe. The closed diamonds with error bars are 
supernova data [^, for which Hq is 67.3 km/s/Mpc based 
on Planck 2013 results [^ . A similar dL-z relation has been 
discussed in entropic cosmology 
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FIG. 2: (Color online). Dependence of the effective equation- 
of-state parameter We on the redshift 2 : for the indicated dis¬ 
sipation rates jl. Zero dissipation corresponds to a nondissi- 
pative ACDM model, whereas ~ 1 corresponds to a fully 
dissipative CCDM model. The background evolution of the 
universe in each case is equivaleiit to that in the fine-tuned 
standard ACDM model because Da = Da = 0.685. The de¬ 
pendence of We on a has been discussed in Ref. [5l[. 

where H/Hq is from Eq. (I28|) . As shown in Fig. [1] it is 
found that cLl for Da = 0.685 agrees with the supernova 
data. (Da = 0.685 used here is approximately equivalent 
to Da = 0.691 for the recent Planck 2015 results [^ . 
The influence of Da is investigated later.) 

To examine the influence of the dissipation rate, fl is set 
to several typical values, 0, 0.05, 0.1, 0.5, and 1.0 in turn. 
The background evolution of the universe in each case 
is equivalent to that in the fine-tuned standard ACDM 
model. As shown in Fig. [TJ the background evolution 
agrees with the observed supernova data because Da = 
Da = 0.685. 

Before studying the radiation temperature, the evo¬ 
lution of the effective equation-of-state parameter We is 
examined. To this end, We for various values of p, is plot¬ 
ted in Fig. [2 where We is calculated from Eqs. (l33t and 
(fTBl) . In this figure, p — 0 corresponds to a nondissipa- 
tive ACDM model, whereas p = 1 corresponds to a fully 
dissipative CCDM model. As shown in Fig. [2 We for 
p — 0 is always equal to 0 because pe = 0. However, We 
decreases with increasing p. It is found that the dissipa¬ 
tion rate p affects We even if the background evolution 
of the universe is not altered. In addition, We for p > 0 
gradually decreases with decreasing z and eventually ap¬ 
proaches — 1. That is. We is not constant when p > 0. 
The varying value of We indicates that the simple T-z 
relation from Eq. (H51) is not suited to describe the radia¬ 
tion temperature in the present model. Accordingly, Eq. 
(I41|) plays an important role. 

Next, the radiation temperature in the modified dis¬ 
sipative model is examined for various values of p. The 
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FIG. 3: (Color online). The radiation temperature-redshift 
relation. The vertical axis is r/(l-|- 2 ). The continuous curves 
represent the modified dissipative model for p — 0, 0.05, 0.1, 
0.5, and 1. The background evolution of the universe for each 
value of p is the same because Da = 0.685. The symbols with 
error bars are observed data points The original 

data represented by the green squares, blue triangles, and red 
circles are from Refs. [JQ, [I], and [a-[Tll]. respectively. The 
recent data represented by the pink diamonds at low redshifts 
are from Ref. [T2|. The data point for 2 = 0 is from Ref. [J]. 
For the observed data points, see Table [I] and Ref. 0 . Note 
that several error bars shown here are asymmetric because the 
corresponding errors shown in Table [T] have been asymmetric. 

radiation temperature T in the late universe is numeri¬ 
cally calculated from Eq. 6ID- For this calculation, the 
redshift is varied from 2 = —0.98 to 4, i.e., the normal¬ 
ized scale factor is varied between a = 50 and a = 0.2. 
To compare with observations, the data points are taken 
from Refs. ilMl. (The values are summarized in Table 
m For details, see, e.g.. Ref. [llj .) In Fig. [2 the vertical 
axis is T/{1 + z) so that a horizontal line can correspond 
to a linear law T = ro(l + z) in equilibrium. From this 
figure, it is found that the radiation temperature T for 
p = 0 obeys the linear law and is consistent with the 
observed data. However, with increasing dissipation /2, 
T gradually deviates from both the linear law and the 
observations. In particular, T for p = 1 is far off the ob¬ 
served data points, due to the nonequilibrium dissipative 
processes. Accordingly, a fully dissipative universe is con¬ 
strained even if the background evolution of the universe 
is the same. On the other hand, T for p = 0.05 agrees 
with the observed data points. This agreement implies 
that a weakly dissipative universe fits the observed T- 
2 relation. The weakly dissipative universe is discussed 
later. 

In the preceding discussion. Da = Da = 0.685. Ac¬ 
cordingly, the background evolution of the universe is 
the same. Consequently, a weakly dissipative model is 
consistent with observations. Finally, to examine the in¬ 
fluence of both p and Da, a likelihood analysis for the 
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TABLE I: The observed temperature-redshift relation. The 
measured values are based on the analysis of a fine structure 
of atomic carbon the Sunyaev-Zel’dovich effect , the 

rotational excitation of CO molecules and a sample 

of X-ray selected clusters [l^. The temperature at a = 0 is 
taken from Ref. [J|. 


z 

T (K) 

Refs. 

0.000 

2.725 ± 0.002 


0.023 

2.72 ±0.10 


0.152 

2.90 ±0.17 

[8] 

0.183 

2.95 ±0.27 

[8] 

0.200 

2.74 ± 0.28 

[8] 

0.202 

3.36 ±0.20 

[8] 

0.216 

3.85 ±0.64 

[8] 

0.232 

3.51 ±0.25 

[8] 

0.252 

3.39 ±0.26 

[8] 

0.282 

3.22 ±0.26 

[8] 

0.291 

4.05 ± 0.66 

[8] 

0.451 

3.97 ±0.19 

[8] 

0.546 

3.69 ±0.37 

[8] 

0.550 

4.59 ± 0.36 

[8] 

0.042 

2.857 ±0.018 

[1^ 

0.077 

2.953 ±0.021 


0.123 

3.072 ± 0.026 


0.169 

3.162 ±0.020 

[12] 

0.222 

3.326 ±0.015 

[12] 

0.274 

3.495 ±0.016 


1.729 


[llj 

1.774 



2.038 


[1^ 

2.418 

9.15 ±0.7 

[9,1^ 

2.690 

10.5l°;« 

[10,1^ 

1.777 

7.2 ±0.8 

m 

1.973 

7.9 ±1.0 

m 

2.337 

10 ±4 

m 

3.025 

12.113,2 

m 


radiation temperature is performed. (A similar analysis 
was performed in Ref. [5ll| . in which a growth rate for 
clustering was examined.) For this purpose, p, and Oa 
are treated as free parameters. The chi-squared function 
for the radiation temperature then becomes 


29 

xrt(^a, m) = 

i=l 


^ohsi^i^ Oa,/ i) 




(47) 


where T'obs(-Zi) and Tcai{zi,Q.f^,p) are the observed and 
calculated radiation temperatures, respectively, and 
is the uncertainty in the observed temperature. The ob¬ 
served data points (numbered f = 1 to 29) are summa¬ 
rized in TablelH For the likelihood analysis, Oa and p are 
sampled in the range 0 to 1 in steps of 0.005. Therefore, 
negative dissipation rates are not considered. Using 



FIG. 4: (Color online). Contours of the normalized likeli¬ 
hood Lrt in the {(^a, p) plane for the radiation temperature. 
The contours of Lrt for 0.2, 0.4, 0.6, and 0.8 are plotted. The 
likelihood function is normalized using the maximum value, 
obtained for (Ha, p) ~ {^a, 0). Note that the maximum value 
is the same for all values of Cl a when p = 0. 


from Eq. gZl), the likelihood function Lrt is [4l| 

Lrt oc exp(-x^T/2)- (48) 

For simplicity, Lrt is normalized. Note that Ua = 0 
and Ua = 1 have not been sampled, in order to avoid a 
division by zero when T is calculated from Eq. (HU- Here 
Ha = 1 corresponds to = 0 because Hm = 1 — Ha. 

FigurelUplots the contours of the normalized likelihood 
Lrt in the (Ha, p) plane. In this figure, (Ha, p) = (Ha, 0) 
corresponds to a pure ACDM model for Ha. Consider 
the contours of Lrt for Ha ~ 0.7. It might be seen that 
low dissipation {0 < p ^ 0.01) corresponds to high Lrt 
regions when Ha ~ 0.7. That implies that a weakly dis¬ 
sipative Universe is compatible with the current data. A 
similar result has been discussed previously (Fij i , in which 
the growth rate for clustering (related to structure for¬ 
mations) was examined. Accordingly, a weakly dissipa¬ 
tive universe is likely proposed. This expectation may be 
confirmed by the bottom panel of Fig. [5] which shows the 
contours of the normalized likelihood L not only for the 
radiation temperature but also for the growth rate. To 
confirm the expectation in more detail, a joint likelihood 
analysis is performed later. 

In the bottom panel of Fig. [Sj the normalized likeli¬ 
hood Lgr for the growth rate is taken from Ref. [Flj . In 
that reference, a likelihood analysis for the growth rate 
in a modified entropic-force model was performed using 
the same method. The modified entropic-force model is 
equivalent to the modified dissipative model considered 
here. Therefore, Lgr for the growth rate in the present 
model is the same as that in Ref. [Fl|. 

In addition, using the same method, the normalized 
likelihood for the distance modulus (related to supernova 




















0.1 
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FIG. 5: (Color online). Top: The normalized likelihood Lsn 
for the supernova as a function of Cl a- Bottom: Contours of 
the normalized likelihood L for both the radiation tempera¬ 
ture and the growth rate. To calculate Lsn, the Union 2.1 
set of 580 type la supernovae is used [^. (For details, see 
the text.) The likelihood Lqr for the growth rate is from Ref. 
[5l| . whereas Lrt for the radiation temperature is replotted 
from Fig. [T] For clarity, the contours of L = 0.2 and 0.8 are 
plotted in the bottom panel. 


data points) is discussed. In the present study, we call 
this the normalized likelihood Lsn for the supernova. For 
the supernova data points, the Union 2.1 set of 580 type 
la supernovae (up to redshift z = 1.414) is used M- The 
chi-squared function is 


580 


Xsn(^a, Ho) — 


i=l 


Uohsi^i) Ucali^i: Hq) 


,.SN 


and the distance modulus ^ is defined as 
fi = 5 log (11 + 25 


(49) 


(50) 


where the luminosity distance is from Eq. (1451) . Keep 
in mind that the distance modulus /x is not the dissipa¬ 
tion rate fl. In this analysis, Cl a and Hq are treated as 
free parameters and sampled in steps of 0.001 and 0.02, 
respectively. That is, the normalized likelihood Lsn does 
not depend on the dissipation rate fl. Consequently, Lsn 
in the {C1a,Ho) plane is obtained. Using Lsi^{CIa, Hq), 
the maximum value at each Cl a can be determined. The 
value is plotted as a function of CIa in the top panel of 
Fig. [SI In this panel, the maximum value is obtained at 
CIa = 0.7221°;°^? (with Hq = 70.04t°;|^ This result is 
consistent with that examined in Ref. |2l| . 

A weakly dissipative universe likely describes the ob¬ 
served data, as mentioned previously. To confirm this 



FIG. 6: (Color online). Contours of the normalized like¬ 
lihood Ltotai for the joint likelihood analysis in the (Ha,/!) 
plane. The likelihood function is normalized using the maxi¬ 
mum value, obtained for (IIa, p) = (0.725, 0). The contours of 
1(T, 2a, and 3a confidence levels (corresponding to —2 In Ltotai 
equal to 2.30, 6.16, and 11.81, respectively [^) are plotted. 
That is, the contours correspond to Ltotai = 3.17 x 10”^, 
4.60 X 10~^, and 2.73 x 10“®, respectively. 


expectation, a joint likelihood analysis is performed, us¬ 
ing the three normalized likelihood functions. For the 
joint likelihood analysis, a combined likelihood function 
Ltotai is defined by 

itotai = Lrt X Agr x Lsn- (51) 

In this analysis, CIa and ft are sampled in steps of 0.005. 
The obtained Ltotai is normalized. The contours of Ltotai, 
corresponding to Icr, 2a, and Str confidence levels, are 
plotted in Fig. [5] The region surrounded by the likeli¬ 
hood contours is close to zero dissipation, i.e., fl = 0, 
when CIa ~ 0.7. In fact, the maximum value of Ltotai is 
obtained for (Ua,/x) = (0.725,0), where an upper limit 
of the dissipation rate is /i ~ 0.01 at the Icr confidence 
level. In this sense, a weakly dissipative model proposed 
here is very similar to the standard ACDM model. 

As has been noted in Refs. [tHzI, adiabatic pho¬ 
ton production does not affect a blackbody spectrum if 
radiation fields are not thermalized. Also, Chluba has 
recently examined the influence of a CMB (cosmic mi¬ 
crowave background) spectral distortion and adiabatic 
photon production processes on the T-z relation in de¬ 
tail j^. Consequently, it is found that the photon pro¬ 
duction process does not affect a blackbody spectrum 
except when the process has a very special energy de¬ 
pendence. In contrast, the CMB spectral distortion af¬ 
fects the CMB temperature. The CMB distortion can 
be constrained using COBE/FIRAS limits. (For details, 
see Ref. [l^.) Accordingly, if adiabatic photon produc¬ 
tion for CCDM models is assumed, a weakly dissipa¬ 
tive model discussed here should be further constrained 
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because the CMB distortion is neglected in this study. 
That is, the dissipation in the universe is expected to be 
smaller than the above mentioned one especially when 
the radiation temperature is discussed. Similarly, the 
recent data 


_ _ Il3|, l23| imply a more weakly dissipa¬ 
tive model. For example, the simple T-z relation of the 
form T = To(l -I- has been investigated in the 

Planck 2015 results [1^. In that reference, one finds 
P = (0.2 ± 1.4) X 10“^, where a recombination red- 
shift oi z = 1100 is adapted. Substituting z = 1100, 
P = (0.2 ± 1.4) X 10-3, and Tq = 2.725 ± 0.002 K 
(from Table HI into the simple T-z relation, we have 
T = 2965 ^ 3028 K. Therefore, when z = 1100, up¬ 
per limits of the dissipation rate can be estimated from 
this result, using m- The upper limit is approxi¬ 
mately jl = 0.018 7^. This constraint is consistent with 
the confidence levels shown in Fig. |6l 


VI. CONCLUSIONS 

The radiation temperature-redshift relation has been 
examined in a dissipative universe. A phenomenologi¬ 
cal modified dissipative model has been developed that 
includes two constant terms, assuming a homogeneous, 
isotropic, and spatially flat universe. The model is equiv¬ 
alent to an extended ACDM model in a dissipative uni¬ 
verse. Therefore, it behaves as if a nonzero cosmological 
constant A and a dissipative process are operative. A 
general radiative temperature law has been used 

in the model, to deduce a T-z relation for a dissipative 
universe. That relation has been computed in the late 
universe as a function of a dissipation rate ranging from 
p = 0, corresponding to a nondissipative ACDM model, 
to p = 1, corresponding to a fully dissipative CCDM 
model. 

The results confirm that the T-z relation for p = 0 
obeys a linear law in equilibrium. However, the calcu¬ 
lated radiation temperature T gradually deviates from 
the linear law with increasing /2, even if the background 
evolution of the universe is not altered. In particular. 


T for p = 1 is nonlinear because the effective equation- 
of-state parameter We varies with time in a dissipative 
universe. In contrast, T for low p agrees with observa¬ 
tions when Da = 0.685 (i.e., the background evolution 
is equivalent to that of a fine-tuned pure ACDM model). 
This agreement indicates that low dissipation describes 
the radiation temperature-redshift relation. That is, the 
dissipation rate is constrained by the observed T-z rela¬ 
tion, even if density perturbations are not treated. The 
present study thus provides new insights into a dissipa¬ 
tive universe. 

To examine the influence of Da, a likelihood analysis 
has been performed. A \ow-p high-DA universe and a 
high-/l 1ow-Da universe have high likelihoods, consistent 
with previous work (Fij i in which a growth rate for clus¬ 
tering (related to structure formations) was examined. 
However, higher and lower Da values are inconsistent 
with supernova data. Accordingly, a weakly dissipative 
universe {p ^ 0.01) for Da ~ 0.7 is a viable scenario. To 
examine this scenario in more detail, a joint likelihood 
analysis has been performed. Consequently, a weakly 
dissipative model proposed here is found to be very sim¬ 
ilar to the standard ACDM model because the expected 
dissipation rate is small. Interestingly, recent works of 
the radiation temperature [mil, [ 23 L ItII imply a more 
weakly dissipative model. However, the properties of a 
weakly dissipative model differ from those of a nondissi¬ 
pative ACDM model. Therefore, further observations are 
necessary to determine whether a low-dissipation model 
is valid. It should be noted that a fully dissipative CCDM 
model (for which p = 1) agrees with observations of the 
growth rate if a negative sound speed [11] and the exis¬ 
tence of clustered matter are assumed. 
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